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Abstract 

■We discuss the novel electronic properties of graphene under an external periodic scalar or vector potential, and the analytical and 
numerical methods used to investigate them. When graphene is subjected to a one-dimensional periodic scalar potential, owing to 
the linear dispersion and the chiral (pseudospin) nature of the electronic states, the group velocity of its carriers is renormalized 
highly anisotropically in such a manner that the velocity is invariant along the periodic direction but is reduced the most along the 
perpendicular direction. Under a periodic scalar potential, new massless Dirac fermions are generated at the supercell Brillouin 
zone boundaries. Also, we show that if the strength of the applied scalar potential is sufficiently strong, new zero-energy modes 
may be generated. With the periodic scalar potential satisfying some special conditions, the energy dispersion near the Dirac point 
becomes quasi one-dimensional. On the other hand, for graphene under a one-dimensional periodic vector potential (resulting in 
a periodic magnetic field perpendicular to the graphene plane), the group velocity is reduced isotropically and monotonically with 
|the strength of the potential. 
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1. Introduction 

Graphene is a single atomic layer of carbon atoms arranged 
in a honeycomb structure. The electronic states in graphene 
obey a unique linear energy dispersion relation near the Fermi 
energy [Fig. 1(a)], and they possess an additional quantum 
number called pseudospin which describes the electron's proba- 
bility amplitudes at the two different sublattices of carbon atoms 
forming graphene JJ]. These behaviors are similar to those of 
massless neutrinos in relativistic quantum physics except that 
Jhe role played by the actual spin of the neutrinos is now re- 
placed by the pseudospin in graphene. 

When graphene is subjected to a slowly varying nanoscale 
external periodic scalar or vector potential, the quasiparti- 
(cles in graphene show even more interesting physics. Re- 
cently, there have been several studies on the electronic pror 
pities _of graphene under either a periodic scalar potential 




under a periodic vector poten- 
or under a periodic comiga- 



itial n 
"tion ll 

Graphene superlattices are not just theorists' dream but have 
been experimentally realized. Superlattice patterns with period- 
icity as small as 5 nm have been imprinted on graphene through 
electron-beam induced deposition of adsorbates [22], triangu- 
lar patterns with ~10 nm lattice period have been observed for 
graphene on metal surfaces 123, 24, 25, 26. 27. 28, 29], and peri- 
odically corrugated graphene sheet has also been reported BOll . 
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Periodically patterned gates can provide another route for mak- 
ing graphene superlattices. 

In this paper, we review the electronic properties of charge 
carriers in graphene under an external periodic scalar or vector 
potential. Especially, we focus on one-dimensional (ID) peri- 
odic potentials [e.g.. Figs. 1(b) or 1(c)] for simplicity. How- 
ever, many of the essential findings discussed here are appli- 
cable to two-dimensional (2D) periodic potentials. We also 
discuss the methodologies used in the analytical and numeri- 
cal calculations. 

The rest of the paper is organized as follows. In Sees. |2]and|3] 
we present the analytical derivation of the energy-momentum 
dispersion relation near the Dirac points (original or newly gen- 
erated) in graphene under an external periodic scalar potential 
and vector potential, respectively. In Sec. |4] we present the 
details of our numerical calculations used in studying the ef- 
fects of stronger perturbing potentials. In Sec. |5] we discuss the 
emerging massless Dirac fermions and quasi- ID modes under 
a strong external periodic scalar potential. Finally, in Sec.|6] we 
summarize our findings. 

2. Graphene under an external periodic scalar potential 

In this section, through analytical calculations, we show that 
when a ID periodic scalar potential is applied to graphene: (i) 
the group velocity of the massless Dirac fermions is anisotrop- 
ically renormalized in momentum space in an unexpected fash- 
ion, and (ii) new massless Dirac fermions are generated at the 
supercell Brillouin zone boundaries ||^. 
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Figure 1: (a) Schematic diagram of graphene. Inset: the linear and isotropic 
energy dispersion near one of the Dirac points in graphene. (b) A ID graphene 
superlattice formed by Kronig-Penney type of scalar potential V{x, v) = V{x) 
periodic along the x direction with spatial period L. The potential is Uo/2 in the 
grey regions and - Uo/2 outside. Inset: energy dispersion of charge carriers in 
this graphene superlattice. The energy dispersion along any line in 2D wavevec- 
tor space going through the Dirac point is linear but with different group veloc- 
ity. For a particle moving parallel to the periodic direction, the group velocity 
is not renonnalized at all whereas that for a particle moving perpendicular to 
the direction of periodicity is reduced the most, (c) A ID graphene superlattice 
formed by Kronig-Penney type of vector potential A(x,y) = Ay{x)y periodic 
along the x direction with spatial peiiod L. The vector potential Ay(x) is Ao/2 
in the grey regions and -Ao/2 outside. Inset: energy dispersion of charge car- 
riers in this graphene superlattice. The group velocity around the Dirac point is 
reduced isotropically . 



We consider a situation where the spatial variation of the ex- 
ternal periodic potential is much slower than the inter-carbon 
distance so that inter-valley scattering between the K and K' 
points in the Brillouin zone may be neglected 13 ll 13211 . We 
shall further limit our discussion to the low-energy electronic 
states of graphene which have wavevector k -H K close to the K 
point, i.e. with |k| <K |K|. 

There are two carbon atoms per unit cell in graphene, form- 
ing two different sublattices. Hence the eigenstate of charge 
carriers in graphene can be represented by a two component 
basis vector The Hamiltonian of the low-energy quasipar- 
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tides in pristine graphene in a pseudospin basis, 
\ .... . / 1\ . / 



/k-r 



and 



e'^'^, where 



and symbolically represent 



1 / ' \ / — \ 1 ; 

Bloch sums of ./r-orbitals with wavevector K on the sublattices 
A and B, respectively, is given by ll33ll 



(1) 



where vq is the band velocity and the cr's are the Pauli matrices. 
The eigenstates and the energy eigenvalues are given by 



V2 



1 



and 



E%) = shvok, 



(2) 



(3) 



respectively, where i = + 1 is the band index and fit is the angle 
between k and the direction. 

We now assume that a ID scalar potential V(x), periodic 
along the x direction with periodicity L, is applied to graphene 
[Fig. 1(b)]. The Hamiltonian H then reads 



H = hvo (- 



■ydy + I V(x)/hV0) , 



(4) 



where / is the 2x2 identity matrix. Next we perform a similarity 
transform, H' - U[HU\, using the unitary matrix 



U^ 



-iaixfjl 



1 

v! 

where a(x) is given by 



ia{x)/2 _g/a(x)/2 



e 



ia{x)/2 



a(x)^2 [ 
Jo 



V(x') dx' /Hvo . 



(5) 



(6) 



Here, without losing generality, we shall assume that an appro- 
priate constant has been subtracted from V(x) and that V(x) has 
been shifted along the x direction so that the averages of both 
V(x) and a{x) are zero. The transformed Hamiltonian H' takes 
the form 



H' = hvo 



-idx 

,-ia(x)g^ 



(7) 



A similar transform has been used to study nanotubes under a 
sinusoidal potential [34, 35]. 

We are interested in the low-energy quasiparticle states 
whose wavevector k s p -h Gm/2 [where Gm = m(2jT/L)x = 



mG X is a reciprocal vector] is such that |p| <K G. In this case, 
we could treat the terms containing dy in Eq. (|7]) as a perturba- 
tion. Also, to a good approximation, H' may be reduced to a 
2x2 matrix using the following two states as basis functions 



1 




^,(p+G„,/2).r ^jjj 



„/(p-G„,/2)-r 



(8) 



Here, we should note that the spinors | q j '^^'^ ( ) '^^^ 

have a different meaning from ^ q j '^'^'^ | ^ 

In order to calculate these matrix elements, we perform a 
Fourier transform of e'"^ 



where the Fourier components /;[y]'s are determined by the 
periodic potential V(x). We should note that in general 



< 1 , 



(10) 



which can directly be deduced from Eq. The physics sim- 
plifies when the external potential V{x) is an even function and 
hence a(x) in Eq. (|6]l is an odd function. If we take the com- 
plex conjugate of Eq. ^ and change x to -x, it is evident that 
//[y]'s are real. General cases other than even potentials are 
discussed in Ref. l'^. For states with wavevector k very close 
to G,„/2, the 2x2 matrix M whose elements are calculated 
from the Hamiltonian H' with the basis given by Eq. ([8]) can be 
written as 

M = hvQ [pxCTz + fmPyO-^^ + hvf, otG/2 ■ / . (11) 

After performing yet another similarity transform M' - 
uImUi with 



1 



1 1 
-1 1 



(12) 



we obtain the final result; 



M' - hvo (pjCTt + /„, /?yO"y) + hvi) m G/2 ■ I . (13) 
The energy eigenvalue of the matrix M' is given by 



^^(P) = Shvo ^pI + l/miVy + m G/2 . 



(14) 



Equation ( fT4l l holds in general and not only for cases where the 
potential V(x) is even i^. The only difference of the energy 
spectrum in Eq. Q from that in Eq. (fT4l i. other than a constant 
energy term, is that the group velocity of quasiparticles moving 
along the y direction has been changed from vo to |/,„|v'(). Thus, 
the electronic states near k = Gm/2 are also those of massless 
Dirac fermions but having a group velocity varying anisotrop- 
ically depending on the propagation direction. The group ve- 
locity along the x direction is unchanged independent of the 
potential. Moreover, the group velocity along the y direction is 
always lower than vq [Eq. ( fT0] i1 regardless of the form or mag- 
nitude of the periodic potential V{x) as schematically depicted 
in Fig. 1(b). 



We have thus shown that other than the original Dirac points, 
new massless Dirac fermions are generated around the super- 
cell Brillouin zone boundaries, i. e. , the case with non-zero m 
values in Eq. (fT4l l. It has also been shown that these newly gen- 
erated massless Dirac points are the only available states in a 
certain energy window if graphene is subjected to a 2D repul- 
sive periodic scalar potential having triangular symmetry ||6|]. 

One more thing to note is that in graphene under an external 
periodic scalar potential, a generalized pseudospin vector can 
be defined and used to describe the scattering properties be- 
tween eigenstates; and especially, back-scattering processes by 
a slowly VMying impurity potential are suppressed as in pristine 
graphene | 



(9) 3. Graphene under a ID external periodic vector potential 



Now we move on to the case where a ID vector potential 
M.x,y) - Ay(x)y is applied to graphene [Fig. 1(c)]. We show 
through a novel transformation relation between scalar and vec- 
tor potentials [18] that, unlike the electrostatic case, the group 
velocity of charge carriers in graphene under a ID periodic vec- 
tor potential is renormalized isotropically . 

The superlattice Hamiltonian, following the Peierls substitu- 
tion, is given by 



H = Svo {^ia-xdx - icTydy - o-yeAy(x)/hcj , 



(15) 



where e is the charge of an electron (e < 0) and c is the speed 
of Ught. The time-dependent Dirac equation then reads 



ih— = hvo (^-icTxdx — icTydy — cryeAy(x)/hcjt// . 



dt 



(16) 



Writing the wavefunction as if/(x,y;t) = e '^'^'^ e''^^'^' (fix), the 
Dirac equation becomes 

E(p(x) = hvQ (^-icTxdx + (Tyky - a-yeAy(x)/hcj ip(x) . (17) 

Now, if we multiply Eq. (fTTl l by dy on both sides, define <f'(x) = 
U3 ip(x) with 



t/3 



1/1 1 



V2\ 1 -1 



(18) 



and use the relations U^, - U'-^ - U^^, UiO-yU^ = -cr,,, and 
U^CTxU-i - <j,, Eq. (fTTl l becomes 

E'if'ix) = hVQ {-iCTxdx + CTy k\, + I V'{x)lhvQ) if'ix) . (19) 
Here, we have defined 

E' = -ihv^)ky , (20) 

k'y^iE/hvo, (21) 

and 

V'(x)^-ie(vo/c)Ay(x). (22) 

Equation (fT9] l is thus equivalent to the Dirac equation with a pe- 
riodic scalar potential in Eq. dU except that now the variables 
are imaginary numbers. Using analytic continuation [Ts^, we 
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obtain the energy-momentum dispersion relation in magnetic 
graphene superlattices from that in electrostatic graphene super- 
lattices. For the states near the original Dirac point, as shown in 
the previous section, the energy dispersion in graphene under a 
periodic scalar potential is given [from Eq. (IT4l 1 by 



(23) 



£,(k) = shvo ^kl + |yi,|2A:2 , 
where, according to Eq. 

/(, = ^ exp {i ^ V(x') dx'^ dx . (24) 

Plugging Eqs. (EB, and ^ into Eq. we obtain 

(25) 



1 



where 



= l j| expjjl ^A,(x')fl'x'Jjx. (26) 



Therefore, from Eq. (125) . we find that the group velocity in 
graphene under an external periodic vector potential (corre- 
sponding to a perpendicular magnetic field) is renormalized 
isotropically in the k^ - ky space l lTlllSll even though the exter- 
nal vector potential profile is highly anisotropic in the x-y plane 
[Fig. 1(c)]. Note from Eq. (|26]l that 



l/o'l > 1 



(27) 



regardless of the form of the vector potential Ay{x), i. e. , the 
group velocity in graphene under an external periodic vector 
potential is always reduced. This result can in fact be used as a 
special case to understand the predictions of velocity reduction 
in metallic carbon nanotubes and gap reduction in semiconduct- 
ing carbon nanotubes under a magnetic field 1.36. .37J . 

4. Numerical calculation 

If one wants to find the energy eigenvalues and eigenfunc- 
tions of graphene under an external periodic potential with 
wavevector k not very close to the supercell Brillouin zone 
boundary centers (Gm/2), one has to resort to numerical cal- 
culations. Such numerical calculations have led us to the dis- 
covery of solutions corresponding to new branches of massless 
Dirac fermions iH that are not found in the analytical calcula- 
tions discussed in the previous sections. 

The scattering amplitudes arising from the periodic scalar 
and vector potentials between eigenstates of pristine graphene 
[Eq. ©I, using the Hamiltonians in Eq. (|4| and in Eq. ( ITsT l, are 
given by 



and 



(^°j/y(r)|<A'j,k,) 

«k|-evo/c-cr,A,(r)|^^,^,) 
evolc ■ Zg k (^'e''^' - ^e-'*") A,,(G) 5k',k-G 



(28) 



(29) 



respectively. Here, G is a superlattice reciprocal lattice vec- 
tor and y(G) and Ay(G) are the corresponding Fourier compo- 
nents of the external periodic scalar and vector potentials, re- 
spectively. Therefore, the energy dispersion and eigenstates of 
the quasiparticles in a graphene superlattice are obtained non- 
perturbatively within the single-particle picture by solving the 
following set of linear equations 

{E-El^ c(5,k) 

= 2.',G 5 + ss' e'^^^'-^^'>) y(G) c(s', k - G) (30) 

for graphene under an external periodic scalar potential, and by 
solving the following set of linear equations 

{E-El^ c(5,k) 

= i:.',G 5 (^'e'"^' - se-''-)Ay(G.) c{s\ k - G) , (31) 

for graphene under an external periodic vector potential, where 
E is the superlattice energy eigenvalue and we have used 
Eqs. (O, (I28]l, and The amplitudes c{s, k) and c{s' , k - G) 
indicate the mixing among diff'erent unperturbed quasiparticle 
states of pristine graphene. 

Note that the scattering amplitude methods presented here 
are applicable to 2D graphene superlattices in general and not 
just to ID periodic systems. 

5. Emerging new massless Dirac fermions in a strong exter- 
nal periodic scalar potential 

If the external periodic scalar potential applied to graphene 
is sufliciently strong, new branches of massless Dirac fermions 
are generated near the original Dirac cone (i. e. , zero-energy 
modes) I'^Sli^. Ul- Note that these new zero-energy modes 



are different from the new massless Dirac fermions discussed in 
Sec.|2] which have higher energies and are generated at the su- 
percell Brillouin zone boundary centers (k = G,„/2) no matter 
how weak the perturbing potential is. As shown in Fig. 2, the 
number of zero-energy Dirac modes increases with the ampli- 
tude of the external periodic scalar potential. These new zero- 
energy modes could have distinguishable signatures in quantum 
Hall [8] or transport measurement [9]. 

The new zero-energy modes are generated when the applied 
periodic scalar potential y(x) has both even and odd symme- 
tries [Fig. 2 or Figs. 3(a) and 3(b)]. When the odd symmetry is 
broken [Fig. 3(c)], new massless Dirac fermions are generated 
but the energy position could be different from zero [Fig. 3(d)]. 
But, if the even symmetry is broken [Fig. 3(e)], new massless 
Dirac fermions are not generated [Fig. 3(f)]. However, even 
in these broken-symmetry cases, the signatures of the modi- 
fied electronic bandstructure could be captured in, e. g. , Landau 
level measurements ^\. 

It is worthwhile to focus on the conditions on the exter- 
nal periodic potential under which the number of zero-modes 
jumps, as shown by the steps in Fig. 2(b) and Fig. 2(e). We 
call these systems at the jump special graphene superlattices 
(SGSs) iSi]. In an SOS, the group velocity along the ky direc- 
tion vanishes. Figure 4 shows that the energy dispersion in an 




Figure 2: (a) Schematic diagram of a Kronig-Penney type of scalar potential applied to graphene given by l/o/2 for < A' < L/2 and - U()/2 for L/2 < x < L with 
lattice period L. (b) Number of Dirac points (not including spin and valley degeneracies) in a graphene superlattice versus Uo. (c) Electron energy in units of 
(= hvo/L; for example, if L = 20 nm, ££ = 33 meV) versus wavevector near the Dirac point for a graphene superlattice formed by the periodic scalar potential 
depicted in (a) with Uo = 6nEt. (d)-(f): Same quantities as in (a)-(c) for a graphene superlattice formed by a sinusoidal scalar potential V{x) = Vq sin(2ff a/L). In 
(f), Vo = 4.0£l was used. 
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Figure 3: (a) Kronig-Penney type of potential V(x) given by Uo/2 for < a < L/2 and -Uo/2 for L/2 < x < L with lattice period L. (b) Electron energy (in units of 
£/. = hvo/L) versus ky with = in a graphene superlattice formed by the periodic potential in (a) with Uo = 6ff£L. (c) and (d): Same quantities as in (a) and (b) 
for a periodic potential V{x) with a perturbation that breaks the odd symmetry. The perturbing potential /W{x) within one unit cell is given by +10 % of the potential 
amplitude {Uo/2) for L/8 < a < 3L/8 and zero otherwise, (e) and (f): Same quantities as in (a) and (b) for a periodic potential ^(a) with a perturbation that breaks 
the even symmetry. The perturbing potential AI'(a) within one unit cell is given by +10 % and -10 % of the potential amplitude (Uo/2) for L/4 < a < L/2 and for 
L/2 < A < 3L/4, respectively, and zero otherwise. 
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Figure 4: (a) Schematic diagram showing the electronic energy dispersion relations and pseudospin vectors (black arrows) in graphene. (b) Calculated electron 
energy bandstructure in graphene. (c) and (d): Same quantity as in (a) and (b), respectively, for the considered Kronig-Penney type of an SGS [as defined in Fig. 1(b) 
with Uo = 4tol where sl = hva/L]. Red and blue arrows in (c) represent the 'right' and the 'left' pseudospin state, respectively. 
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SGS is quasi- ID and the pseudospin in an SGS is either paral- 
lel or antiparallel to the +kx (or periodic) direction [5]. Because 
of this quasi-lD electron energy bandstructure in an SGS, the 
group velocity of electronic states are almost the same over a 
wide region in momentum space; and this has led to the pre- 
diction that SGSs can be used for electron beam supercollima- 
tion lH. 

6. Summary 

The electronic structure of graphene under a general exter- 
nal periodic scalar potential is modified from that of graphene 
in several highly unexpected ways: (i) the group velocity is 
anisotropically renormalized in momentum space and (ii) new 
massless Dirac fermions are generated at the supercell Bril- 
louin zone boundary. Moreover, when a strong ID periodic 
scalar potential is applied, new zero-energy modes emerge. Un- 
der certain conditions, the group velocity of charge carriers 
along the direction perpendicular to the ID periodic direction 
of the superlattice potential vanishes. In this special class of ID 
graphene superlattices, the electron energy bandstructure and 
pseudospin structure are quasi- ID and these properties can be 
used in collimating the electron flow. With a ID external peri- 
odic vector potential applied to graphene, on the other hand, the 
group velocity of charge carriers near the original Dirac point 
is reduced isotropically. The analytical and numerical methods 
discussed in this paper can be used in further investigating the 
novel properties of quasiparticles in ID and 2D graphene su- 
perlattices. 
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